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We propose to integrate the electro-optic (EO) tuning function into on-chip domain engineered lithium 
niobate (LN) waveguide. Due to the versatility of LN, both the spontaneously parametric down conversion 
(SPDC) and EO interaction could be realized simultaneously. Photon pairs are generated through SPDC, 
and the formation of entangled state is modulated by EO processes. An EO tunable polarization-entangled 
photon state is proposed. Orthogonally-polarized and parallel-polarized entanglements of photon pairs are 
instantly switchable by tuning the applied field. The characteristics of the source are theoretically 
investigated showing adjustable bandwidths and high entanglement degrees. Moreover, other kinds of 
reconfigurable entanglement are also achievable based on suitable domain-design. We believe tailoring 
entanglement based on domain engineering is a very promising solution for next generation 
function-integrated quantum circuits. 

Entanglement is a crucial physical resource for quantum information science and technology. In practical 
investigations and applications, entanglement systems based on photons are widely used owing to the extra 
long decoherence time. To perform various specific tasks, different kinds of entanglements are required, 
such as multi-particle entanglement, nonmaximally entanglement, and mixed state entanglement'"^. For the 
generation of entangled photons, spontaneous parametric down conversion (SPDC) in nonlinear crystals is 
one of the most powerful tools". In SPDC process, a single photon is transformed into a photon pair. However, in 
previous works' it is usually necessary to design complicated optical paths to transform the initially generated 
states into target entangled states. Numerous large-scale bulk optical elements have to be fixed onto sizeable 
optical benches in the laboratory, such as lenses, wave plates and filters'-"''". The physical size and the inherent 
instability of these components hinder the development of more complex schemes, and bring significant chal- 
lenges for the integration of quantum circuits, which needs to be improved in practical applications" '^. 

In this article, we propose a different strategy for on-chip tailoring entanglement through domain engineering 
in a lithium niobate (LN) waveguide. LN is a typical ferroelectric nonlinear optical material, which is simulta- 
neously available for SPDC'"" and electro-optic (EO) effect'''"'. Due to its quite large effective nonlinear coef- 
ficient, LN could serve as entangled photon-pair source of high quality""*. What's more, EO modulation is an 
effective way to manipulate entangled photons. Previous investigations have demonstrated the control of phase" 
or polarization state'''. An integrated quantum relay operator has been established based on EO effect of single 
domain LN crystals with some special electrode design''\ In fact, at the quantum level, EO could be regarded as a 
two-photon interaction process, in which the quantum states of photons could exchange. In this work, we derive 
relevant theories to describe this process. Combining the theories for EO and SPDC, entangled states of desired 
formations could be designed flexibly. Through suitable domain engineering, which is widely used in quasi phase 
matching (QPM) nonlinear optics", the EO and SPDC processes could be effectively combined to generate the 
needed entangled states. Moreover, LN is also a good platform to realize versatile functions including efficient 
infrared photon detection^", negative permittivity^', EO quantum logic gate^^ and lensless ghost imaging^'. AU of 
them could be integrated together toward future practical large scale quantum circuit integration. 

Following, as an illustration of our entanglement design approach, we establish an EO tunable polarization- 
entangled photon pair source in a domain engineered LN waveguide with suitable artificial structures. When a 
suitable voltage is applied, the source produces a pair of entangled photons contains the same polarization, with 
either o- or e- polarization state. The notations o- and e- represent ordinary and extraordinary light in uniaxial 
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anisotropy crystal, respectively. In contrast, if the voltage is turned 
off, the entangled photons bear orthogonal polarization states, i.e., 
one photon is o-polarized while another one has e-polarization. Since 
we may control the polarization state intentionally, the formations of 
entangled states thus could be regulated accordingly, which possess 
considerable potential in modulation of multi-photon entangled 
states. 

Results 

Tailoring entanglement in a domain engineered LN waveguide. In 

recent works^''"^'', integrated quantum circuits established on silicon 
waveguides are adequately investigated. Utilizing the femtosecond 
laser waveguide writing technology, on-chip devices to support and 
manuscript polarization-encoded qubits have been fabricated"'^'^". A 
controUed-NOT gate has been experimental demonstrated with 
switchable entanglement. It is no doubt that silicon waveguides 
have excellent capabilities to confine and propagate light. 
However, limited by the inherent properties of silicon, they are not 
very convenient for fast manipulating and modulating entangled 
photons; such as the modulation scheme based on thermo-optical 
effect^^, the response time still has much room to improve in practical 
applications. In addition, to realize multifunctional integrated 
quantum circuits, entanglements of complicated formations are 
normally required such as the nonmaximally entanglement''^ and 
mixed state entanglement*''"''. Under these circumstances, the LN- 
based system is considered as an alternative choice. Besides excellent 
nonlinear optical properties, EO effect in LN material provides an 
effective way to change the quantum state of photons through the 
photon interaction process. For a LN crystal, if we set its symmetry 
axis as the z-axis and apply a voltage along the y-axis, a variation is 
led into the dielectric constant. It could be written as £ = £(0) + /fe 
with AEjk = — eo'Xsi-EflH^Mg {jk = 23, 32)"". In the equation, y^i is the 
effective EO coefficient, and represents the applied electric 
field. The refractive indices of ordinary and extraordinary light are 
marked as and n^. This additional portion of permittivity 
corresponds to an equivalent polarization item, which is expressed 



as Pi=-'^ e,ay^iEanln]h.jkEk with A^/t = 



1, jk = (23, 32) 
0, else 



From 



d rP-E, we obtain the interaction Hamiltonian for EO 



process as 



0 

-L 



— i{o)' — (a)f\ 



with 



Ffo — 



dydzU7)U^- 



(2) 



In the equation, the quantized expression of electric field is utilized. 
Feo is the overlap integral of waveguide modes expressed as i9(r). The 
function /fxj represents the structures of the inverted domains. The 
EO coefficient is periodically modulated, thus the o- and e-beams are 
coupled due to the periodic index ellipsoid deformation. The phase 
matching condition should be satisfied between the ordinary and 
extraordinary photons to ensure high conversion efficiency. The 
first item of the equation represents the annihilation of an e- 
polarized photon and the creation of an o-polarized photon, while 
the second corresponds to the inverse process. Accompanied by the 
two photon interaction process, the old photon state is substituted by 
new state which is desired to form target entangled state. Choosing 



proper EO processes to modulate entangled photon states supplies 
an innovative approach for entanglement architectures. 

As a simple illustration, the polarization qubit of photon is con- 
sidered. Assuming Hi is acted on a photon with polarization state | o) , 
it would be transformed into polarization state |e) with a certain 
possibility. Referring to the corresponding classical situation"', we 
could obtain a superposition of | o) and | e) through modulating the 
applied voltage. The ratio of these two portions is determined by the 
value of applied voltage. The overall effect could be equivalent to H^o 
= R{!X.)Hi, where J? is a standard rotation operator, and a is the 
equivalent rotation angle. For a biphoton state, the modulation pro- 
cess could be written as 

|*P/)=H£01®H£02|*P,) (3) 

with the matrix expression of the operator as 

sin 0 sin (J) — sin 0 cos (j> — cos 0 sin (f) cos 0 cos ^ 

- sin 0 cos (j> — sin 0 sin (f) cos 0 cos (f) cos 0 sin (fi 

- cos 0 sin (j> cos 0 cos (fi — sin 0 sin ^ sin 0 cos ^ 
cos 0 cos cos 0 sin ^ sin 0 cos ^ sin 0 sin (f> 

where 0 and (p are the corresponding equivalent rotation angles to 
determine the relative ratio of the two photon vector states, which are 
controlled by the applied voltage. It works similar to a half waveplate 
(HWP). However, there are still apparent differences. Assuming the 
input light is a linearly polarized light, a circular polarized beam is 
obtained if the wave plate thickness is only half In contrast, the light 
passing through a PPLN with half-length is still linear-polarized, 
while the orientation angle is also half if the applied field is the same. 

For universal multiphoton entangled state, they could not be gen- 
erated through SPDC directly, so the postprocessing modulations are 
quite critical. In multiphoton entanglement situation, only certain 
formations may lead to potential applications. We need to choose 
proper photons from the initial state F{jyi); 1^2) !^n)}, where 
{|4'„)} represents the photon pairs from SPDC, and the function 
-Fd^')) corresponds to their initial combination formation. After 
suitable modulations are applied based on the effective ffgo. the state 
thus could be transformed into a valuable formation such as GHZ 
state or cluster state. This process could be expressed as 



1*/)= ® <i^{|^i);l*2)- 

|*/c,>{|*„)} 



I^m)}, 



(5) 



where (x) refers to the direct product of all -f/^Q. 

After the corresponding scheme is decided, the ferroelectric 
domains of the LN waveguide could be well-designed to satisfy the 
phase matching conditions of the related SPDC and EO processes. 
Through analyses of the set of [Ak^, the matching vectors could be 
obtained as {2n/A,„}. Accordingly, the corresponding Fourier coeffi- 
cients and the Fourier bases could be derivate. Based on these, the 
original function is recovered through 



(6) 



where /frj represents the actual domain arrangement of the entang- 
led photon-pair source. It is worth mentioning that the correspond- 
ence between {Ak„} and {2n/A^} does not have to be one-to-one. The 
choice of {A„} is flexible, which is an important part of domain 
engineering. 

EO tunable polarization-entangled photon state. As a detailed 
illustration, we consider the generation of an EO tunable 
polarization-entangled photon state with switchable characteristic. 
We consider the processes of photon interactions in a z-cut, x- 
propagating titanium in-diffused LN waveguide. There are well 
defined transverse electric (TE) and transverse magnetic (TM) 
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propagation modes in the waveguide at near-infrared frequencies^'. 
In this work, only the fundamental modes of the TE and TM modes 
are discussed. For convenience, they are marked as o-polarized and 
e-polarized, respectively™. Three photon interaction processes are 
included. They are two SPDC processes, which are expressed as 
Os + e„ Op ^ gj + 0, ip, s, and i represent the pump, signal, 
and idler waves, o and e correspond to ordinary and extraordinary 
polarization); and an EO process takes place for the signal wave, 
namely ^ {e^ ^ oj. If there is no voltage applied, only the 
SPDC processes happen in the waveguide. The power of the pump 
light is launched into the downconversion photons in corresponding 
propagation modes of the waveguide, the generated entangled state is 
type-II-like with the formation |os,e,) + |es,Oi). In contrast, if there is 
an applied voltage, the EO effect is activated. The EO polarization 
rotation and SPDC processes happen simultaneously and are 
coupled everywhere in the crystal. In this case, the down-converted 
photons are equally generated everywhere and affected by the locally 
EO perturbed refractive index ellipsoid. The final obtained 
polarization state, including, e^. Op and o^ are the superposition 
result of all "sub-sources" inside the crystal. Therefore the rotation of 
polarization is realized through coherent superposition of all kinds of 
photons randomly collected in the sample. Corresponding to the 
coupling wave theory, that means the power of the fundamental 
TE mode (Oj) and the TM mode (fij) at the signal frequency 
exchanges with each other. For a given sample length, the 
efficiency is dependent on the interaction strength, which is 
proportionate to the applied voltage. We could always find a 
proper voltage to make the efficiency reach —100%, then the 
entangled state is switched to be type-I-like with the formation 
\Ss,Si} + \os,Oi). Thus the entangled state is switchable through 
controlling the applied voltage on and off 

To satisfy phase matching conditions of the three interaction pro- 
cesses, the domain structures need to be well designed. General 
speaking, we may design three sets of domain periods so that their 
reciprocal vectors are able to compensate these three wave vector 
mismatches. However, this straightforward scheme normally 
induces very thin domains and complicated structures. To solve this 
problem, several approaches have been proposed. For example, for 
the optimized wavelengths and domain structures, the needed 
domain periods could be reduced to two and they are of integer ratio, 
so that the domain structure could be easily fabricated. The schem- 
atic is shown in Fig. 1. Correspondingly, the sign of nonlinear coef- 
ficient x*^' and EO coefficient y change periodically. The modulation 
functions are expressed as d(x) = d3ifi(x)f2(x) and y(x) = 
y5ifi(x)f2(xf, where /(xj = sign[cos(2n/A)] = 'L^Gmexp(ik^x). In 
the equations, d is the substitute of nonlinear coefficient x'-^'' with the 
relationship d = x*^'/2. The effective component of d and y in our 
situation are d^j and y^i, respectively. Thus the reciprocal vectors are 
given as = 2mn/A, with their corresponding Fourier coefficients 
as Gm = (2/mn)sin(mn/2). As a consequence, the expansion formula- 
tion of the coefficients are obtained as"" 

m,n 

(7) 

yW=y5i^G^,ne'^'"'"' 

with 

4 mn ,nn^ 

Gm,„= rsin( — )sin( — ) 

mnn'^ 2 2 

8 

2mn 2nn 

where Ai and A2 represents the corresponding two modulation per- 
iods. Based on these analyses, the equations of QPM conditions could 
be expressed as 




Figure 1 | Schematic of a dual-PPLN waveguide. The Kght and dark blue 
portions represent the positive and negative domains of the PPLN, 
respectively; while the purple portion is the core of the waveguide. The two 
golden strips correspond to the electrodes. A 1 and A 2 are the two periods of 
the structure. Moreover, the polarization directions of o- and e-polarized 
photons are corresponding to y- and z-directions, respectively. 

~ Pp.0 ~ Ps,o~ Pi,e ~^mi,ni 
Nh=Pp,o-lis,e-ko=Km,,n,- (9) 

In these equations, „ (j = p, s, i; cr = o, e) refers to the propagation 
constant of the corresponding waveguide mode, which is obtained 
based on the Hermite-Gauss formulations^'. zl/Jj and Afi2 corre- 
spond to two SPDC processes, while Aji^ corresponds to the EO 
interaction. Kmj,nj (i = 1, 2, 3) are the required reciprocal vectors 
to compensate the phase mismatches. 

Through detailed calculations, we find an appropriate set of solu- 
tions for Eq. (9). The width and depth of the waveguide core are both 
set at 10 |im, and the maximum of index difference is set at 0.003. 
The pump, signal and idler wavelengths are 0.7335 |im, 1.6568 |im 
and 1.3162 |tm, respectively. As room temperature operation is 
always more appreciated for future quantum circuits, the simulation 
temperature is chosen at 25°C. The photorefractive effect might 
bring some influences. However, the photorefractive damage of 
LN is not a simple process. It depends on many factors such as laser 
pulse width, operation wavelength and surface quality. Through 
recent technologies^'^, the damage threshold has been increased 
to an acceptable degree for many applications. Moreover, the MgO 
doping is also an effective way to further avoid crystal damage'*'''". 
In this situation, the corresponding values for the (m, n) series 
are {(m^, «,) = (3, 1), (m^ «2) = 0, -1), (m^, n^) = (1, 1)}. 
The dual-modulation periods are obtained at A^ = 25.84 |im 
and A2 = 154.96 pm, respectively. The ratio of these two periods 
is A2/A1 = 6 and the duty cycle is 0.5, as is shown in Fig. 1. (More 
details about the selecting process could be referred to the Methods 
section). 

To give a quantum mechanics description of our tunable entang- 
led photon-pair source, we derive the state vector of entangled 
photons through the effective Hamiltonian. The total interaction 
Hamiltonian consists of two parts, which is expressed as Hj = 
HspDC + Heo- Hspoc corresponds to the SPDC processes, while 
H^o refers to the EO interaction. Both of them arise from polariza- 
tion, namely, Pf^ji and Peo- In detailed treatments, the pump field 
is usually treated as an undepleted classical wave. The signal and 
idler fields are quantized and represented by field operators. Thus 
the formulations of the pump, signal and idler fields could be written 
as 
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doh 



doji 



ho). 



2£nnlN, 



(10) 



2eoniNi 



where cr represents the polarization state o or e, and nj„ (j = s, i) is the 
corresponding refractive index. and Nj^ are normalization para- 
meters. Impair), i9siT(r) and i9,^(r) correspond to the transverse mode 
profiles of pump, signal and idler field, respectively. 

For SPDC processes, if using the rotating wave approximation, 
HspDC is derived as 



d^rPNi:E = t:n 



rf;,g.-(/V-fto-A. + Jf»,.)x 



dmdwF I — fltgt ,(o.p-».-a,i)t (11) 



-L J) V niniN,eN,o 



+ H.C.] 



with 



foe = Jd:t'£iz,9po(r),9,„(iO<9i^(r) 
d;'ciz.9po(r)5s,(r),9io(r) 



(12) 



For the EO process, the interaction HamUtonian is derivate as above. 
H^Q is expressed as 



Heo = — ysi-Eo Gm,« [ 



(13) 



V ^-^^'"^ 



with 



Feo = dydz3so{r)^%e{i) 



(14) 



The entangled state vector could be obtained through 



|y(0>=exp[(-//K) 



dt(HsPDC+HEo(Ea))]\0). Detailed treat- 



ments thus could be divided into two cases: 

i) when there is no applied voltage, we expand the evolution oper- 
ator to the first-order perturbation, so the state vector is 
expressed as 



|>F) = |1--| dtHspDcimo) 



fee 
— rC' 

: I dte-'<°''-'°'-°''"4aJj|0) 



In the equation, we have 



and 



dxe"^'** = /j(lA/;) with the formation h(x) = ( 



-ix/2 



sinc(x/2). Moreover, the actual SPDC process doesn't merely 
arise at the perfect phase-matching frequencies flj and Qj. As 
the natural bandwidth is v, we set cOj = il^ + v and cOj = — v, 
then the state vector is simplified as 



W = |o)- 



+ 



dvP,MLAp,Jalal\0) 
dvPMLNheoWj,M 



(16) 



with 



Poe — i^LEpodiiGi^iFo, 



Pea — i^LEpodji Gi^ _ iFj, 



(17) 



In the equation, the phase mismatch terms could be expressed as 
^Pooc = Pp.0 - Ps.o - P i.c and A = Pp_„ - P,_, - pi_o-As |v|< 
i3s, Qj, we have substituted wXoJj) by QsiSii) in Poe and P^o- 
ii) For the second case, the voltage is applied. The evolution 
operator should be expanded to the second-order term. The 
formulation of the state vector could be expressed as 

|4') = |l+(-^f| rffzj <i(,T(Hjo(f2)Hsi.Dc(ti))||0) 

" " (18) 



+- 



X [G3, MLA!S„,)P„ aL-J, + G3. -MLAI1„JP, [- ^ 4°Ll}l°) 

In the equation, T is the time-ordering operator. The notation 
zf^oj corresponds to the mismatch of /^^ „ and P^ ^. The functions 
h(LAPooe) and h(LAPoeo) are spatial integration for SPDC, while 
h(L/ipce) and h(LApgJ correspond to the EO processes. There 
are four product terms of effective operators, including 
alusealal, aj,as„aj„fl],, aj„fls,aj„a], and fl>,oaJ,fl;„. The first 
two correspond to actual processes. If they act on the vacuum 
state, the results show flj„aseaj^fl,'„|0) = (oJ^jOso)) 
(fl,,fljj0,,)) (flL|0,„) j = 1 1,„,0,„1,„) and oJ^fl,„flJ„oJJO) = 

ial\Ose)){asoalK))(al\Oie))=\he,Oso,he}- For the last two 



terms, they don't affect the final entangled state because as 
«se|0)=0 and flsojO)=0. Therefore the parallel-polarization 
entangled state is written as 



l*) = |0) + 



+ 



dv 
dv 



dv'P„MLApJh{LAp,Jalal\0) 



dv'PMLAPJhiLAp„Jal/jO) 



(19) 
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with 



0 J m 2 
Poo= — — "3iy5l"s£po£aGl, 1G3, - ipEoFeoL' 



2 V nlNlN„N,„ 



0 J r\ 2 
-Pi;e= asiyjitisEpo^aGl, 1G3, ipEoFoeL 



(20) 



Similarly, we also set cOj = fJj + v and cu,- = f3, — v, while cUj (to,) 
is substituted by (f2,) in Poo and P„. Summarizing the two 
situations above, we rewrite the entangled state vector in a terser 
formation: 



l^> = |0) + ( 



d{Pa)dEa)[ 



dvPoMLAPooeWsoal\0) 



dvP,oh{LAI]oJalal\0}]+{l 



8{EMPa) 



(21) 



dv 



+ 



dv'PoohiLAIUKLAPoJal^Uo) 
dv'P^MLAUhiLAPooeWAm 



In the equation, S(EJ represents the Dirac-(5 function. When 

|-0 + 

= 0, we have 

' ~ ' ~ " \e)s\o) 

fO + 

5(Ea)dEa) should be 1, the 



corresponds to Coe|o)Je), + C, 
nonzero value, the value of (1 - 



S(E„)dEa = 1. The generated entangled state 
If Ea is an appropriate 



0- 



entangled state vector is C„o\o)^\o)- + Cee\e)^\e)j. To ensure the 
EO transformation of the downconverted photons, several cal- 
culations have been done. The crystal length is set at 3 cm, and 
the average pump power is set at 1 W. Based on the coupled 
wave equations for EO interaction"" and those for downconver- 
sion process, we obtain that when is 3.8 X 10^ V/m (Usually, 
the distance between the electrodes used for the LN waveguide is 
10 ~ 1000 |im, the corresponding applied voltage is about 1 ~ 
100 V)'^, the transformation efficiencies are PJiP^e + Pso) ~ 
0.9927 (for \o) p^\o) ^\e) .^\e) ^\e) ,) andPJiP,, + PJ = 0.9932 
(for |o)p— >|e)s|o)i— >|o)s|o),)) respectively. As for the idler 
photon, the polarization state is hardly rotated. That's because 
the PPLN's EO effect is quite special with wavelength selectivity. 
Only lights at the designed discrete wavelengths are affected. The 
corresponding wavelength bandwidth also could be pre- 
designed. For a direct verification, the transformation efficiency 
of polarization of the idler photon is calculated to be about 1%. 
Therefore the corresponding influence could be neglected. The 
orthogonally polarized entangled pair could be transformed into 
parallel polarized entangled pair with a near 100% efficiency. 
Therefore EO tunable entangled states are realized. Further- 
more, due to the intrinsic properties of EO modulation, the 
entangled state is fast switchable. That could be realized in both 
CW and pulsed pump regimes. In the CW pump regime, we 
could just set the voltage at a certain value until different output 
photon states are requested. For a 3-cm long sample, the switch- 
ing time is limited by the propagation time from the entry to the 
exit surface, which is about 0.2 ns. Being different from the CW 
pump, the pulsed regime may give some more constrains. We 
should avoid changing the voltage when the pumped pulses 
propagate right in the sample. Since the highest repetition 
rate of commercial lasers are only around 100 MHz, the time 
gap between cascading pulses should be long enough for EO 
sv«tching. 



Characterization of the tunable entangled photon-pair source. To 

ensure the quality of our tunable entangled photon-pair source, we 
investigate the corresponding properties. Firstly, the spectrum 
character is discussed, which is represented by the modulus 
squares of h-functions. For type-II-like polarization entangled 
photons, the corresponding expressions are \h{LA[iooe)\^ and 
\h{LA[iogo)\^' whUe those for the type-I-like polarization entangled 
photons are simultaneously affected by SPDC process and EO 
interaction, so the formulations are \h{LAp^^)h(LAP^^o)\^ and 
\h(LAflo^)h{LApooe)\^ > respectively. Calculation results are 
presented in Fig. 2. The natural spectral bandwidth of entangled 
photon-pair source is mainly influenced by the group velocities of 
the signal and idler wavelengths. In our situation, these two 
wavelengths are close with each other, so the differences between 



the bandwidths for the two SPDC processes, /. e., Op 



Os + ei. Op 



+ o, are relatively small. From Fig. 2(a) and Fig. 2(b), when the 
waveguide length L is 3 cm, these two quantities are 0.35 nm and 
0.27 nm, respectively. Compared with the previous report™, 
these values are relatively smaller, which is beneficial to the 
improvement of entanglement degree. On the other hand, the 
natural bandwidth increases with the decrease of the waveguide 
length L. For instance, we could see that the value of bandwidth 
increase to —0.9 nm for Op ^ + e, process when the crystal 
length reduces to 1 cm. For Coo|o)Jo),- + Cee|e)j|e),- pairs, as the EO 
process is integrated together, the natural bandwidths are further 
restricted through the corresponding h-functions. If the waveguide 
length L is 3 cm, the bandwidths for ordinary and extraordinary 
photon pairs are 0.20 nm and 0.23 nm, respectively. The 
difference between these two values is smaller, and the spectrum is 
almost symmetrical. 

The anticorrelation dip could be calculated based on the spectrum 
function, which is usually different for type-I and type-II SPDC""". It is 
the Fourier transform of the SPDC spectrum. The corresponding 

formulation is _Rc(t) ~ 1 — (l/2;r) |ft(v)|^ cos (vt)c?v. A narrower 

dip corresponds to a broader spectrum, so the width of anticorrela- 
tion dip is mainly determined by the natural bandwidth of the SPDC 
spectrum. For the orthogonally polarized entangled pair, |/i(v)|'^ 
corresponding to sinc^{LA[iooe/2) and sinc^{LApoeo/2)- We expand 
the phase mismatching to the first nonzero order of v. The 

corresponding formulations are Apooe~ ~iPso~ f^ie)^ ^^'^ N^oeo~ 
~iPse~ Pio^^'- Pja 0 = p, s, i; C) = o, e) represents the dispersion 
parameter. For the paraUel polarized entangled pair, as its spectrum 
is further confined by EO process, the spectrum property is naturally 
different from those of the traditional cases. The |/2(v)|^ functions are 
written as sinc^{LAIiooel'^)'sinc''{LAPoJ2) and sinc^{LA[ioeo/ 
2)-sinc\LAPJ2). Similarly, we have A/J„,«(/4 -/4)v and 
APeo~(l^se ~ l^so)^'- iUustration, the anticorrelation dips corres- 
ponding to sinc^{LAfSooe/2) and sinc^{LAPooe/2)' sinc^{LA[SoJ2) are 
plotted in Fig. 3. 

Based on the discussions above, we calculate the von Neumann 
entropy of our source, in which the influences of Pks (k. 8 = o, e) 
are included. It is defined as S = — tr(ps„(,log2Psuf,)''. Psub represents 
the reduced density operator for the subsystems. For a product 
state, the quantity vanishes. If the state is maximally entangled, the 
value of the entropy is 1. To describe tunable entangled photon-pair 
source, the density operator of the entangled state is expressed as 
follows 



Ps, = 
with 



0 + 



8(Ea)dEa ]\(Po){(Po 



<5(£Jd£j|^p)(<^p|(22) 
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Figure 2 | Normalized output signal spectra corresponding to the two type of entangled states, respectively, (a) Op 



, + e,-; (b) o„^ e, + o,; (c) o„ 



+ e,+ e;, (d) o^- 
L is 3 cm. 



■ fij + O5 + Oi. The corresponding natural bandwidths are 0.27 nm, 0.35 nm, 0.20 nm and 0.23 nm, when the waveguide length 



ko) = 



+ 



dvP„,fi(LA/?„„,)aJ„flI|0) + 
dv f dv'P,MLA[Uh{LAPoeo}alal\0) 
dv dv'P,MLAUKLAPooe>Ul\0) 



dvP,MLAPoeo)alal\0) 



(23) 



The reduced density operator could be obtained through 
Ps = — fnOsi)- Therefore the quantity S is derived as 



•0+ 

0- 
•0+ 



5{Ea)dEa 



5{Ea)dEa 



log2 



0p '"^2 y 0p 



00 



©O 

Qee 



1 / ®™ 



(24) 



log2 



0. 



where ©„.= 



dvPoMLAPooe) 

2 



0p '"^2 y 0^ 



dvp,„h{m„j 



, 0a 



dv'P^h{LApjh(LAp^) 



and 0a, = 



dv 



dv'PMLApjh(LAp^) 



Besides, 0o = 0oe + ©eo and 0p = &oa + ©ee- 

From the equation above, we could see that the value of S is mainly 
determined by 0ij (i, j = o, e). In an ideal situation, i.e., Q^c = ©eo 
and ©00 = 0ee> S would reach its maximal value 1, which corre- 
sponds to the maximal entanglement. If the natural bandwidth is 



narrow, and the ratio Poc/Peo. Poo/Pce is close to 1, ©oe and ©„□ thus 
are almost equal to ©eo and ©ee, which makes S ~ 1 and a very high 
degree of entanglement. On the contrary, wider bandwidth would 
affect the entanglement accordingly. As is discussed above, the nat- 
ural bandwidth in our case is quite small. Following the treatments in 
Ref 38, we simplify the calculation of S into the perfect phase match- 
ing situation. Thus S has the same formulation for two types of 
entangled states. This result verifies that EO process does not affect 




Delay i (nm) 



Figure 3 | The anticorrelation dips Rc(t) are calculated based on 
sinc^(Ljp„ae/2) and sinc^(Ljp„o^2) x sinc^(LJp„^2), respectively. The 
blue solid line corresponds to the former, while the red dashed line 
represents the latter one. The crystal length is set at 3 cm. 
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the entanglement degree. Besides, we could see that is mainly 
influenced by geometry and material parameters from Eq. (21). We 
find that the source has high tolerance for the geometrical variations. 
Until the waveguide modes of and a, are nearly cutoff, S stiU stays 
above 0.99. The degree of entanglement maintains in a high level. 

Discussion 

What we have presented above is merely a simplest case for the two- 
photon entangled state, where the entangled states are generated and 
manipulated through EO modulation. Moreover, this type of system 
also could be used to generate multiphoton entanglement. A straight- 
forward proposal is considered based on the switchable two-photon 
entanglement source proposed above. For our device, if the voltage is 
turned on, the generated entangled state would be |ee) -|-| oo). We 
may increase the pump power to improve the simultaneous genera- 
tion possibility of double photon pairs in SPDC processes. In this 
situation, we may obtain four states with the same possibility, 
namely, |ee)|ee), |oo)|oo), |ee)|oo) and joo)|ee). To demonstrate 
the entanglement, we'd better erase influences of |ee)|oo) and 
|oo)jee). That could be realized through a post-selected detection 
scheme, which is widely used in energy-time entanglement. After 
the generated photon pairs pass through a single domain LN wave- 
guide, the arrival time of o-polarized and e-polarized photons would 
be different due to the birefringence. If the frequencies of the down- 
converted photons are close and chosen around telecom band, the 
frequency dispersion would be small enough to be neglected''. 
Therefore |ee)|ee) and |oo)|oo) could be detected the same time, 
while there would be time difference in detection of photons in 
|ee)joo) and |oo)|ee) states. With a narrow time window that is 
usually used in energy-time entanglement*, a four-photon GHZ 
state could be selected and detected. Moreover, with the help of 
on-chip optical elements, such as, mode analyzer, directional coup- 
ler, Mach-Zehnder interferometer and Y-coupler"'^^'''^, entangled 
state of more than four photons and preselected entanglement 
scheme are designable. 

Moreover, with the help of artificially designed domain structures, 
the nonmaximally entangled state^'^ could also be easily realized, 
which can find applications in the loophole-free tests of Bell inequal- 
ities*' and quantum metrology with the presence of photon loss''"'. 
The LN waveguide could be divided into several functional portions 
along the propagation direction"'^^. We may set the first portion to 
regulate the amplitude ratio of o-polarized and e-polarized pump 
light through EO interaction. Since both the o-polarized and e-polar- 
ized pump lights are involved, two corresponding SPDC processes, 
namely, Op^ Os + O; and ep ^ e^ + C; may happen simultaneously in 
the second portion. They are dominated by the nonlinear coefficients 
A22 and d33, respectively. In this case, the QPM conditions also could 
be satisfied for these two processes through suitable domain design. 
The non-maximally entangled state thus could be obtained that is 
expressed as |os,o,) -\-e\es,ei). The relative ratio 8 of amplitude maybe 
controlled by the voltage applied on the first portion. As a con- 
sequence, an EO tunable non-maximally entanglement source is 
realized. The mixed state entanglement is not hard to be generated 
based on nonmaximal entanglement. A usual procedure is to create 
an entangled state with desired degree of entanglement, modulate the 
density matrix and then introduce decoherence""'. The first step 
could be realized through non- maximal entanglement with a suitable 
entanglement degree, and then we could utilize the EO polarization 
rotation effect to modulate the non-maximally entangled state to 
generate off-diagonal terms in the corresponding density matrix. 
To introduce decoherence into the system, experiments using bulk 
optical elements have been demonstrated based on the spacial prop- 
erties of the SPDC emission cone'' '. However, in waveguide systems, 
it's not convenient to use spatial decoherence directly. We could 
consider the temporal decoherence. Through introducing birefrin- 
gent retardation, the vector | oe) and \ eo) could be distinguished from 



1 00) and \ee) via the relative arrival time difference between two 
photons. Tracing over time, coherence between distinguishable 
terms in the density matrix could be erased, thus specific off-diagonal 
terms become smaller to form the target mixed entangled state. 

It's worth mentioning that besides the uniform domain periods we 
used all above, the domain structure of LN or similar materials could 
be elaborately prepared in more complicated patterns, such as quasi- 
periodic structures, which could provide complex reciprocal vectors 
to compensate the multi-process phase mismatching with high effi- 
ciency"". We believe the domain engineered nonlinear waveguide is 
really a very promising platform for quantum integration circuits, 
which deserves more in-depth studies in both theories and 
experiments'". 

As for practical applications of our device, the quantum cryp- 
tography protocol may be a potential area, which is the most real- 
izable quantum information application nowadays''^. Based on the 
fast switching entangled state, our device provides a different strategy 
for quantum cryptography. Instead of switching the complementary 
analysis basis at the users' (Alice and Bob) places, we could switch 
between the states directly at the source while keeping the analysis 
basis always oriented along the same direction. Although this strat- 
egy may cause certain constrains, an advantage is that the operation 
conveniences for users may be improved. Another possible applica- 
tion example would be the classical and quantum communication 
without a shared reference frame'"'. Alice may encode a logical qubit 
as |0i)=(|HH) + jVy))/v/2, \li)={\HV) + \VH))/s/2, and use it 
to communicate with Bob. If one applies a rotation to both of the 
photons, \Qi) is unchanged, i.e., it is invariant under arbitrary rota- 
tions of the reference frame around the propagation axis. Although 
j li) wUlbe changed, it is stiU orthogonal to \ Qi). Therefore without a 
shared reference frame or the pre-alignment of Alice and Bob's 
polarization basis, a projective measurement on jO^^) by Bob will 
distinguish the two logical bases. Furthermore, with small postpro- 
cessing modifications based on on-chip optical elements mentioned 
above"'^^'"", our source can switch between the generation of two 
entangled states with different symmetries (\HH) + \ W))/\/2 and 
( I HV) — I VH) ) / \/2, both of which are invariant under the rotation of 
the reference frame"". Thus a symmetry measurement by Bob will 
distinguish between the two states. The similar idea can also be 
applied to the error correction protocol in quantum computations. 
Moreover, in most of the integrated platforms to date, the quantum 
information is encoded in the spatial modes^"*"^'. Manipulating the 
polarization of photons on chip allows us to encode and decode 
classical and quantum information in the polarization degree of 



0.01 




Signal wavelength (urn) 



Figure 4 | The values of J(X) against the signal wavelengths when 'Lp is 
equal to 775 nm. The blue solid line corresponds to ( ni}, = (1,1), while 
the red dash line corresponds to ( m^, H3) = ( 1 , — 1 ) . Only the blue line has 
crossover point with the zero line, so ( 1 , 1 ) is possible solution for the 
equation A(a) = Ap^ik) — Gj_±j = 0. 
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Signal Wavelength ().im) Signal Wavelength ().im) 

Figure 5 | In the situation where the subscript of Jp^ corresponds to (inj, 113) = (1, 1), The function is tuned by the pump wavelength and the 
operation temperature, respectively. The green dashed line marks J(l) = 0. (a) Values of zlf^j change with the signal wavelengths and the pump 
wavelengths, (b) Values of A(?.) changes with the signal wavelengths and the operation temperature. From the — 0 lines, suitable pairs of pump, 
signal and idler wavelengths could be chosen for specific situations. 



freedom. Fast switching means the shorter operation time of qubit, 
which is a desirable feature in practical applications. We believe over 
1 GHz modulation bandwidth would be easily achieved in this case. 

In summary, we theoretically propose an effective approach to 
tailoring entangled photons through the combination of SPDC and 
EO effects. Based on the domain-engineered LN waveguide, an EO 
tunable polarization entangled photon-pair source is realized. It 
could selectively generate orthogonal- polarization or parallel-polar- 
ization photon by adjusting the applied voltage. The bandwidths and 
entanglement degrees of both types of entangled states are of excel- 
lent performances. In comparison with the bulk or in-fiber beam 
splitting, polarization handling, and wave plate elements, we pro- 
posed an integrated solution possibly on a single monolithic LN chip. 
This on-chip integration would be very desired in future quantum 
information related applications since it could be more compact, low 
loss and stable. In addition, the EO tunability is another attractive 
feature that supplies extremely fast manipulation on entangled 
photons, which is superior to the fixed wave plate solution. 
Moreover, the manipulation of some other entangled states are also 
discussed based on our proposals, including the multi-photon 
entangled state, the mixed entangled state and the nonmaximally 
entangled state, which show some attractive advantages. 
Introducing EO polarization modulation into domain engineered 
nonlinear waveguides offers opportunities for on-chip manipulation 
of quantum states and fast reconfigurability of the dedicated 
quantum chips, which is very desired for practical applications in 
future quantum circuits. 

Methods 

To find a proper set of pump, signal and idler wavelengths, they are scanned carefully. 
Firstly, when the pump wavelength is fixed, calculation results reveal that A (i j and 
^1^2 (for SPDC processes) are close but much greater than Aj^^i (for EO interaction). 
Thus they are assumed to correspond to the solutions {3, ±1) for {m, n), respectively. 
For certain pairs of signal and idler frequencies, the values of Ai(a) and A2( a) could be 
calculated through linear superposition of Aj^iik) and A§2QX a^iid then we compare 
the linear combination (1, ±1) of 2k/ Ai and 2k/ A2 with A^3 to see if they are 
matched. The corresponding expression is A(X) — Aji^iX) — Gi ±}. For the wave- 
length satisfying A (a) — 0, three phase mismatches could be compensated 
simultaneously. 

In this study, we choose the pump wavelength at 775 nm. The values of A (a) 
against the signal wavelengths are plotted in Fig. 4 (A.p — 775 nm). It is seen from the 
figure that only the linear combination {1, 1) of 2k/ A j and 27i:/yl2 provides a solution 
for the equation A(?.) — 0. Thus Gjj G3 -1, Gij are selected as the compensations of 
A^j, A§2 a^iid Aji^, respectively. Moreover, the phase matching conditions could be 
adjusted by the pump wavelength and the operation temperature, which are corre- 
spondingly shown in Fig. 5(a) and Fig. 5(b). To ensure the ratio between A j and is 
an integer, we finally choose the pump, signal and idler wavelengths at 0.7335 |j.m, 
1.6568 )j.m and 1.3162 |j,m, respectively. In this situation, A2 ~ 6A1, where the 
operating temperature is set at 25^C. For quantum communication applications, the 



operation wavelengths might be preferred to be chosen from the telecom-band. That 
could be easily realized by finding the solutions of Eq. (9) accordingly, as is shown in 
Fig. 5(a). It is clearly from the figure that the whole telecom C band (1.53 [xm- 
1.57 |im) and L band (1.57 |im-1.61 \xm) could be selected for operation. 
Furthermore, the operating temperature could also be conveniently adjusted 
according the simulation results as shown in Fig. 5(b). 
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